DIFERENCIAS E INTERPOLACIONES 


PRINCIPIOS JENERALES 

El objeto de una tabla matematica es el de permitir en- 
contrar el valor de una funcion correspondiente a cuales- 
quier valor del argumento variable. Siendo imposible tabu¬ 
lar la funcion para todos los valores del argumento, tenemos 
que construir la tabla para ciertos valores especiales sola- 
mente. Per ejemplo, en las tablas de los senos i cosenos los 
valores de las funciones son dados para valores del argu¬ 
mento que difiere uno de otro en un minuto. 

El procedimiento de hallar los valores de las funciones 
correspondientes a valores del argumento intermediary en- 
tre los que son dados se llama «interpolation 

En Astronomia con mucha frecuencia se hace uso de las 
interpolaciones, debiendo por tanto los astronomos familia- 
rizarse con estos calculos. 

Debemos hacer notar en primer lugar que, en verdad, 
ningun procedimiento de interpolacion puede ser aplicable 
a todos los casos que se presenten. Por el hecho de que: 

al numero 2 corresponda el logaritmo 0,30103 
i al numero 3 « « 0,47712 

no estamos justificados para sacar cuaiesquiera conclusion 
respecto a los logaritmos de los numeros comprendidos entre 
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2 i 3. Por lo tanto ? una o mas hipotesis sera necesario siem- 
pre adoptar como base de cualesquier sistema de interpola* 
cion. Las hipotesis que se adoptan siempre son estas dos: 

1. Que } suponiendo que el argwmento varie uniformemente 7 
la funcion variard de acuerdo con alguna lei regular. 

2. Que esta lei debe ser estudiada en los valores de las fun¬ 
done# dados por la tabla. 

Estas hipotesis se aplican al procedi mien to para form al¬ 
ias «diferencias », cuya manera de hacerlo se vera en el 
ejemplo que sigue, aplicado a los logaritmos senos de los 
arcos comprendidos entre 20° i 30°: 
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La columna f nos da la diferencia entre cada dos valores 
consecutivos de las funciones. Estas diferencias se llaman 
l. a diferencias. 

La columna f” nos da la diferencia entre cada dos prime- 
ras diferencias consecutivas. Estas se llaman #. a diferencias . 


Del mismo modo, los numeros de las columnas que siguen, 
una vez escritos, se llaman 3. as , 4diferencias , etc.. .. 

Ahora, continuando el orden sucesivo de las diferencias 
veremos que ellas se hacen cada vez mas pequenas o que 
converjen hacia cero, esto nos indica que los valores de las 
funciones siguen una lei regular i por lo tanto, la l. a hipo- 
tesis que hemos hecho es entonces aplicable. 

A fin de poder aplicar la interpolacion debemos suponer 
entonces que los valores intermedios de la funcion siguen la 
misma lei. La verdad de esta presuncion debe ser confirma- 
da de algun modo antes que podamos interpolar con exac- 
titud materaatica, pero en la practica debemos suponerla 
xacta a falta de una razon que nos demuestre lo contrario. 

F.FECTOS DE LOS ERRORES EN LOS VALORES DE LAS FUNCIONES 

En la Tabla A se puede notar que si continuamos el orden 
de las diferencias mas alia de la 4. a empezaran a auraentar 
i a ser cada vez mas irregulares. Esto proviene de las im 
perfecciones de los logaritmos a causa de la omision de las 
decimales mas alia de la o. a cifra. 

A1 notarse que las diferencias ernpiezan a ser irregulares, 
podamos juzgar si esta irregularidad proviene de errores en 
los numeros orijinales, los que deben ser correjidos, o de 
pequenos errores que provienen necesariamente de la falta 
de decimales en los logaritmos. 

La gran ventaja de las diferencias es que cualesquier 
error, aunque sea pequeno, en las cantidades diferenciadas, 
a menos que siga una lei regular, saltara a la vista por las 
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diferencias. Para demostrar la razon de esto, investigare- 
mos mas adelante que efecto producen loserroresen lasfun- 
eiones dadas sobre los ordenes sucesivos de las diferencias. 

Estudiemos ahora el caso en que las funciones sigan una 
lei regular. 

Sea, por ejemplo, x = a 3 , es decir que: x = f (a). Demos 
a a los valores sucesivos de 0 a 10, elevemos en seguida es- 
tos numeros a la 3. a potencia i construyamos la tabla de las 
diferencias: 


tabla B. 





Diferencias 



Fun cion 
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fiv 


a=0 

a 3 =0 

“l~i 




a=l 

a s = l 

+7 
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+ 6 


a = 2 

QO 

II 

10 

cS 

+19 

+ 12 

+ 6 

0 

a=3 

a 3 =27 

+37 

+ 18 

+ 6 

0 

a=4 

a 3 =64 

+61 

+ 24 

+ 6 

0 

a=5 

a 3 =125 

491 

+ 30 

4“ 6 

0 

a=6 

a 3 =216 

+ 36 


0 



+127 


+ 6 


a=7 

a 3 =343 

+169 

+ 42 

+ 6 

0 

a=8 

a 3 =512 

+217 

+48 

+ 6 

0 

a=9 

a 3 =729 

+271 

+ 54 



a=10 

a 3 =1000 
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En la Tabla B las 3. as diferencias son constantes i las 4.^ 
diferencias se reducen a 0 (cero), luego las funciones son 
exactas. 

Estudiemos ahora con detencion las tablas de las diferen¬ 
cias para deducir en seguida las formulas de interpo 
lacion. 

Copiemos nuevamente la l. a parte de la tabla A. 

tabla A 





Diferencias 


A rg amen to 

Fun cion 





f 

f” 


/ IV 



Log. sen 20° 

9.53405 

+ 2028 




21° 

9.55433 

+ 1925 

— 103 

+ 8 


22° 

9.57358 

+ 1830 

—95 

+ 8 

0 

23° 

9.59188 

+ 1743 

—87 


0 

24« 

9.60931 
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—4 






+ 2 






—1 


Un atento examen de estas diferencias nos permitira en- 
contrar para cualesquier argumento su valor, espresado en 
terminos de las diferencias. Asi, porejemplo, el log sen de 
23°=log sen 22° + /’, lo que se puede tambien espresar de 
este otro modo: 
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/ (a + nw)=/ (a) + /' 

Basandonos en esto, insertamos a continuacion una tabla 
que espresa, de una manera jeneral, las relaciones que guar- 
dan entre si los terminos de los ordenes sucesivos de las 
diferencias correspondientes a una funcion que varie si- 
guiendo una lei regular. 
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Facilmente se coraprendera la construccion de esta tabla. 
La suma de los terminos unidos por las lineas diagonales 
producen los terminos indicados por las lineas verticales. 
Asimismo, es facil comprobar la exactitud de dichas espre- 
siones en el ejemplo numerico dela Tabla B. En dicho ejem- 
plo se ve que partiendo de las 3. & s diferencias (f ”), que son 
constantes, se puede encontrar el valor de un ndmero n de 
un orden cualquiera de las funciones. 

Supongamos aliora que en las funciones orijinales haya- 
mos cometido un pequeno error i para discufir el caso en la 
forma mas jeneral que se pueda presentar, imajinemonos 
que hayamos cometido un error de 0.5 en uno de los mime- 
ros ? estando los otros correctamente, pues no es posible su- 
poner que varios mimeros esten equivocados. Construyendo 
la tabla de las diferencias obtendremos: 




Diferencias 


Funcion 








/“ 

/ m 

fiv . 

f Y 

0 

0 

0 

+0.5 

+0.5 

—2.5 

0 

+0.5 

+0.5 

—1.5 

-2.0 

+5.0 

0.5 

—0.5 
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+0.5 
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0 


0 
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En este caso, el valor maximo de la diferencia del orden 
n es: 1.5 en las 3.as diferencias; 3.0 en las de la 4. a ; 5.0 en 
las de la 5. a , etc. 
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Su espresion jeneral es: 

1 n (n—1) (n—2).(n— s+1) 

2 1.2.3...s ~~ 

en que n es el or den de las difereneias i 

n n—1 

S — T ° ~ 2 ~’ 

segun si n es par o impar. Asi: 



1 

2 


i siendo 

2 i / 11 1 n (n—1).(n—s + 1) 

S ” 2 — 1 . J " 1.s 


j_ v 2(1).(2-1+1) 1_ _2 

2 A 1.... 1 ~ 2 A 1 


i siendo 

„ 3-1 , n (n—l)(n—2).. (n—s + 1) 

2 " J 2X 1.2. s 



3 (2) (l).... (3-1+1) 

1.2.1 



i siendo 

4 „ /iv 1 w n(n— l)(n—2)(n—3). . (n— s + 1) 

S ~ T _ y — “1.2.3.s 

i_v 4 (3) (2 ) ( 1). • ■ (4-2 + 1) 1 4.3 

~ 2 A 1.2.3.2 2 2 

i siendo 
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J_ n (n—1) (ii-2) (n—3) (n—4). . (n—s + 1) 

2 A 1.2.3.4.s 

- J_v & W(3)(2)(1). .(5-2+1 )_Lv^i -5 

2 S' 1.2.3.4. 2 2^2 ° 

Este es el caso jeneral, pero en lapractica las difereneias 
pueden resultar dos o tres veces mayores sin que esto im- 
plique forzosamente un error mayor de 0.5 en los numeros 
escritos. 

Tenemos entonces la siguiente regia jeneral para juzgar 
si una serie de numeros siguen verdaderamente una lei uni¬ 
forme: 

«Diferenciemos las series hast a alcanzar un or den de dife- 
« rencias en que los signos + i — se alternen o sigan uno a 
« otro irregularmente . 

«Si ninguna de las difereneias de ese orden, espresada en 
« unidades de la ultima cifra de los decimoles 9 excede del 
« Umite. 

n (n—1) . (n—s+1) 

1.2.3.. s ~ 

«< — esto es, del valor del mayor coeficiente hinomio del orden 
« n, — los numeros dados se consideran que siguen una lei 
« regular, i por lo tanto que son exactos dentro de una unidad 
« de la ultima cifra. 

«Si algunas difereneias exceden de este Umite, su cuociente 
« multiplicado por el coeficiente hinomio de mas arriba, debe 
« considerarse que indica el error mdximo de que probable - 
« mente estd afectado el numero opuesto.» 

De esa manera podremos, con gran seguridad, encontrar 
un error aislado en una serie de numeros. Supongamos, por 
ejemplo, un error de dos unidades en alguno de los numeros 
de la serie que sigue: 
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Construyendo la tabla de las diferencias encontraremos 
que las 4 raayores diferencias del orden 5.° seran: — 10, 
-l- 20, —20,+ 10, lo que nos permitira encontrar al momento 
el numero erroneo i juzgar de la magnitud de su error. Un 
error cerca del principio o fin de la serie de los numeros, 
cuyas diferencias han sido calculadas, no puede ser descu- 
bierto por las diferencias a menos que sea mui conside¬ 
rable. 

Si, por ejemplo, el primero o ultimo numero esta equivo- 
cado en una unidad, el error en cada orden de las diferen' 
cias sera solo 1, como facilmente se puede ver en el ejemplo 
que sigue: 
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DIFERENCIAS 

FUNCION 

Z 1 

f u 

fill 
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— l 




0 

0 


— 1 

etc.., 

0 


0 




Es solamente en aquellas diferencias que estan en o cer- 
ca de la misma linea que los numeros crecen en la forma 
que hemos indicado mas arriba. Pero al principio i fin de 
las series no podemos determinar esas desigualdades. 

Examinando las diversas tablasde diferencias, que hemos 
indicado, vemos que n numeros tienen n—1 primeras dife¬ 
rencias, n—2 segundas diferencias i asi sucesivamente, dis- 
minuyendo el mimero en 1 unidad con cada orden sucesivo. 
Por ejemplo, una serie de 7 numeros tendra 6 primeras di¬ 
ferencias, 5 segundas diferencias, 4 terceras diferencias, 
etc., etc* Por lo tanto, si el ntimero de f unciones dadas no 
excede del indice que espresa el orden de diferencias que se 
deban formar, no se podra obtener ninguna conclusion. 

Lo que se ha dicho aqui sobre la exactitud de los nume¬ 
ros, cuando las diferencias siguen una lei regular, debe en- 
tenderse como aplicable a errores aislados solamente. Si to- 
dos los numeros estuvieran sometidos aun error quesiguiese 
una lei regular, este error no seria descubierto por las di¬ 
ferencias, porque, por la naturaleza misma del caso, las lil- 
timas solo indican desviaciones de alguna lei regular. 
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FORMULAS FUND AMENT ALES DE INTERPOLACiON 
Fdrmula de Newton 

Elijamos una serie de ntimeroseuyas diversas diferencias 
aparezcan como en el cuadro que sigue, el cual contiene una 
seleccion de series que preceden i siguen a una de ellas: 
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En esta tabla se puede ver que todas las diferencias, cu- 
yos terminos son equivalentes, se -encuentran colocadas so- 
bre una misma linea horizontal. Esto sirve para indicar en 
que linea cae una diferencia de cualquier orden. 

Asi, todos los terminos (a + 2) estan en una misma linea 

5 

horizontal i los terminos (a -t- —) se encuentran media linea 
mas abajo. 

Las diferencias que siguen dan una idea clara de la for¬ 
ma como se construyen todas las diferencias de la tabla 
anterior: 


/(a-2 


i t 





i 

/(a + (n + 1) w - /(a + nw)=/ 


(a + (n + 



f (a + n + 


w 


■f (a + ( n — 4-lw = 


f (a + nw )=/ ( a + 1) 


Trazando en la Tabla D, a partir de la funcion orijinal, 
una linea diagonal descendente, obtendreraos todas las dife¬ 
rencias que se emplean en la formula de iuterpolacion de 
Newton, cuya espresion jeneral es la siguiente: 


I 

/ (a + nw) =/(a.) + n f 



n (li — 1) 
2 


II 

/ 


(a + 1) -j- 


n(n — 1 (n —2) ™ 

" 1.2.3 ”f 


a +Y) + 
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n(n — 1) (n — 2) (n— 3) ^ , 

— 1 - ~y* 2'"3 4 -" / IV (a + 2) + .... 


En algunos testos alemanes los coeficientes binomios: 

n ( n — 1) n (n— 1) (n — 2) 

1.2 ; 1.2.3 

se representan de esta otra manera: 


t) > (t 


Entonces, la formula de Newton se puede espresar tarn* 
bien de este modo: 

/(a + nw) = /(a) + n/‘ ( a + If) + ( a + 1 ) + 


\±, 

1 5 ! 


l/ m a + 


j / IV (a + 2) + 


(*)/( 


& + 2 1 ^ 


Aplicaciones de la formula de interpolacion de Newton 

Supongamos que se hayan calculado las coordenadas de 
un cometa para las fechas que se espresan a continuacion: 
Julio 23, Julio 27, Julio 31 i Agosto 4 i deseamos conocer sus 
coordenadas exactas para Julio 24. 

Construyamos la Tabla de las diferencias: 
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Fijaudonos en esta tabla, en la Tabla D de las diferen- 
cias i aplicando la formula de Newton obtendremos los si- 
guientes valores: 
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A esta formula se le puede dar una senciila trasformacion sacando a n como factor comun en los 
diversos terminos del 2.° miembro, i entonces la formula se trasformara en esta otra: 
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Escribamos nuevamente la formula de Newton trasforma - 
da, solo hasta el 3. er termino, aplicandole el valor de n: 

y[ix + 0’,4w) =y^a)-f0’,4 a + — ^ 

O’4—1 f 11 0’,4— 2 f nI / 3 \ 

+ — L 2 - 1 f ( a + 1)+ 3 \f ( a + ”) + •••• 

-=r/(a) + 0’,4 |f ^a + 

f n f ni/ 3 N 

-0’,3{y (a+l)-0,533|y (» + —)+.... 

Ahora, apliquemos a esta ultima formula los terminos de 
las diferencias comprendidos en la linea diagonal descen- 
dente de la tabla de arriba i obtendremos: 


= 9,53405 + O’,4 f 2028—0,3 / 103—0.533 x 8 


824 


1 


—4.27 


log sen 20°24’ =9.54229 


+ 32. 


.18 { 


—107.27 


j 2060,18 
824,07 ( 


«De esta manera, por medio de la formula de interpola- 
« cion de Newton, se pueden encontrar los logaritmos inter- 
« mediarios de aquellos que se dan en las tablas i este es el 
« procedimiento empleado para construir las Tablas de lo- 
« garitmos que todos conocemos i de las cuales hacemos 
« tanto uso en nuestros calculos». 

Formula de interpolacion de Gauss 

Mirando con atencion la Tabla D i relacionandola con al- 
guna de las tablas de diferencias numericas quehemos dado, 
veremos que las diferencias que usa Newton en su formula 
tienen las siguientes equivalencias: 
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ii ii hi/ 1 \ 

/ (»+l)=/ (»)+/ ( a +y) 


hi / 

? I a + 

y ' 



+ 



(a) 


pv(a+2) = fiv(a4 1) + fv ^ j 

a + 4-)+ 

Con un ejemplo numerico cualquiera, es mui facil com- 
probar que estas formulas son enteramente exactas. 

A las diferencias que resultan, corresponderan asimismo 
los siguientes coefieientes binomios: 



i+-|)=f iv ( a )-f fv ( a +4-) 

fv.(a+l)=fiv(a)+2 fv( a -f -|-j 


Para f’ 

( a +4) 

mrrespondera.n 

n. n—1 


Id. f” 

(a) 

■a. 2 


Id. f 


I i( i n(n-l) 

n(n—1) (n—1) 

I ld - 1.2 + 

1.2.3 


i sacando factor coraun a 

n (n —1) 


despues de haber multiplicado por 3 el numerador i deno- 
minador de la primera fraction, obtendremos: 


Para f | a + —jcorrespondera. ... 

n (n—1) ( ol „ | (n + 1) n (n —1 
' 1.2.3 i 3 + n “ 2 I = -E2"3- 
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Id. f iv(a) correspondera 11 ^ 11 ^ ——-f 

1. Z . o 


n(n—l)(n—2)(n—3) 

1 . 2 . 3.4 


i sacando factor comun a 

n (n—1) fn 2) 
1 . 2 . 3.4 


despues de haber raultiplicado por 4 el nunaerador i denomi- 
nador de la priraera fraccion, obtendremos: 

n (n—1) (n—2) ( „ . _ „ ) _ (n+1) n (n—l)(n—2) 

1.2.3.4 1 I - 1.2.3.4 

Luego la formula que resulta, agrupando las diferencias 
con sus respectivos coeficientes binomios, conocida con el 
nombre de «Formula de Gauss» ; se representa de esta raa- 
nera: 

f (a+nw)=! (a) + nf | a + -~j+ f”(a) 

+ (11+ l 1) 2 n . 3~ 1) ( a + y) + W 

(n + 2) (n+l)n (n—1)(n-2) / J_\ 

^ 1 . 2 . 3 . 4.5 \ ^ 2 ) 

Reproduciendo la tabla D veremos que las diferencias 
empleadas por Gauss, siguen una linea zig zag que baja i 
sube: 











Tabla D 
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Tomemos ahora de esta misma tabla las diferencias que 
se eneuentran en una misma linea horizontal i eoloquemos 
al lado de ellas sus diferencias equivalentes. 

Entonces obtendremos: 


f 'f a + T) = f (“-¥)+'” (a) 

f'”(a+i)-r a-+f'v (a) 

A estas nuevas diferencias corresponderan tambien los 
siguientes coeficientes binomios: 


A f ( a — 

\ 

* f” (a) 
»f’”| a — 
» f IV (a) 



corresponded.. .. n 


n(n —1) n (n+1) 

n + — 1 —J . = —(2 + n—1)—n ±-Z-L 


(n + 1) n (n—1) 
1.2.3 


(n-(-2) (n + 1) n (n—1) 
1.2.3.4 


Ahora, agrupando estas diferencias con sus respectivos 
coeficientes binomios obtendremos la 2. a Formula de Gauss, 
cuya espresion jeneral es la que sigue: 

f (a + nw)=f (a)tnf! a-( a ) 

)+ (n + 2>( 1 n 2 + 3 1) ° (n ~‘ ) 


(n+l)n(n-l) _ 1 

TTT 1 ^ a - 


(n + 2) (n + 1) n ( n —1) (n-2) ^ / 1\ 

1 . 2 . 3 . 4.5 1 [ 2 ) 

Trasladandonos de nuevo a la Tabla D veremos que 
Gauss en su segunda formula solo usa diferencias alterna- 
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das que suben i bajan, a partir de la funcion orijinal, i que, 
para liacerlas resaltar mas a la vista, las hemos unido con 
una raya. 

Hagamos ahora una recopilacion de las tres formulas que 
ya hemos determinado: 

Formula cle Neivton 


f(a+nw)=f (a) + n f’ j a +-i-j — f” (a+l) 

■ n (n—1) (n—2) ( 3) 

+ 1.2.3 * l 1" 2 j 

+ ( n ~ 3 . L f iY (a+ 2)+.... 


Formulas de Gauss 


(D I 


1 f(a-fnw) — f (a)-f-nf |a + ~j -f- i ^ ( a ) 


(n-f 1) n (n— 1 ) f „. 

1.2.3 


( a+ t) 


+ 


( 2 ) 


f (a-f mv)=f (aj+n f ^ a — A ) 


(n +1) n 
1.2 


f" (a) 


, (a+l) n (n—1) 1 \ . 

\ +-TO-' l a -“J + . 

La formula de Newton solo usa diferencias situadas en 
una lmea diagonal descendente, como puede verse en la Ta- 
bla D. 

La primera formula de Gauss solo contiene diferencias si¬ 
tuadas en una linea zig-zag que baja i sube. 

La segunda formula de Gauss solo contiene diferencias 
situadas en una linea zig-zag que sube i baja. 

Para interpolar hacia atras se usa siempre la segunda for¬ 
mula de Gauss, eambiando n por — n. Veamos ahora como 
queda dicha formula despues de hacer esta modifieacion. 
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Formula de Gauss jpara la interpolacion en sentido inverso 

(n — 1) n 


f (a — nw)=f (a) — n f ’ | a-- j + 

i f »(„ __L1 - (n+1) n (n — 
l 2 J 1.2.3. 


(n—l)n(n -f 1) f „ , 
1.2.3 


f” (a). 
1) (n—2) 


f!v(a) 


El ejemplo que sigue, da una idea mas clara de la manera 
como se trasforma el coeficiente de la tercera diferencia i 
ayuda tambien a comprender como se trasforman los otros 
coeficientes binomios 


(n+1)—n(—n— 1) 
1.2.3 


(n—1) n (n+1) 
1.2.3 


Notese una particularidad de esta formula que es: la de 

tener los signos-1- alternados en los terminos del 2.° 

miembro, lo que permite recordarla con facilidad. 

Vamos ahora a dar una formula jeneral que sirva al mis* 
mo tiempo para la interpolacion en sentido directo como en 
sentido inverso. La formula es esta: 


yja±» wj= y ( a ) ± ny(a ± i)+ a J 


+ 


(n—1) n (n-f-l) 


1.2.3 


J 


( a± 3-) 


+ 


(n—2) (n—1) n (n +1) 
24 



El signo superior sirve para la interpolacion en sentido 
directo i el signo inferior para la interpolacion en sentido 
inverso. 

Esta ultima formula puede tambien representarse de la 
siguiente manera: 


f (a ± n w) — f (a)±n 


{ 



(n—1) 
2 



(n+1) 


3 
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|f’” (a ± jliv (a) ± | fv (a± l)±. 

La tablaque sigue nos hace ver que las diferencias usa- 
das en las dos formulas de Gauss abarcan una faja circuns- 
crita por dos lineas horizontales trazadas arriba i abajo de 
lafuncion orijinal, quedando en la mediania de esta faja las 
diferencias pares de a . 


FUNCION 

DIFERENCIAS 

f 

f” 

f ))7 1 

fiv 

f (a—2 w) 

«'<»-!) 




f (a—w) 


f” (a—1) 




f(a-i) 


*”(•-?) 


f (a) 

'> + i) 

f” (a) 

_ 

r'(a+-l) 

fiv (a) 

f (a+w) 


f ” (a+1) 



f (a-f-2 w) 
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EJERCICIO 

Calcular el logaritmo seno de 22° 24’ con la ultima fdrmula 
de interpolation de Gauss, interpolando en sentido directo 

Formula: 

f (a+nw)=I (a)+n j!’ (a j(”(aH—Si- 1 

{f”>+4)+T^ {'"(a) 

Como la interpolation, que vamos a hacer, es en sentido 
directo solo hemos tornado el signo +. 

Construyamos ahora la tabla de las diferencias, contenien- 
do solo aquellas que se necesitan en la formula. 

TABLA DE LAS DIFERENCXAS 


ARGTO. FUNCION 


sen 21° 


sen 22° 


sen 23° 


9.57358 


DIFERENCIAS 


(+1925) 


+ 1830 


—95 


(-87) 


(+ 8 ) 


+ 8 


fiy 


0 
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Nota. —Las diferencias que no entran en la formula estan 


colocadas entre parentesis. 

Repitiendo la formula i aplicandole a sus diversos termi- 
nos las diferencias numericas de la Tabla que antecede i sus 
valores a los respectivos coeflcientes binomios, obtendre- 
mos: 





n + l 

3 



24’ 

siendo n= — — = 0/4 


60 
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o 

X 


CM 

I 

O 


o o 

* X 




o o 




TOMO CXXYII 


56 


log sen 22° = 9,57358 
log sen 22o 24' = 9.58101 
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Formula de Stirling 

Eepitamos la tabla jeneral de las diferencias: 


Fnnciones 

Diferencias 

r 

r 

V” 

f (a— w) 

! 

/ f (a) 

( a_ 4) 

t” (a —1) 

, f” (a) 

ri ■ - t) 

1 ! f («) 

1 f’” (a) 

r (“+■?) 

’”(■ + t) 

1 

1 

)i 1 (“ + i)_ 

f (a + w) 

V { IV 

v ) 

f” (a + 1) 


En esta tabla, las diferencias encerradas dentro de un pe- 
queno cuadro: 

-|~J, r (a), f” | a + (a) 

represen tan los promedios de las diferencias entrelas cuales 
estan intercaladas. 

Eepitamos ahora las dos formulas de Gauss: 
f (a+nw)=f (a) + n f [ a + y] + ” ( a ) 
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+ flLt l 2 »(a = i> r (.+ t] . —^5,- ' I IV (a> 


(n+l)nfn— l)(n—2) 
24 ” r 


f (a-fnw)=f (a)-f-n f'| a-^ 


r (a) 


(n+l)n(n—1) 
6 


( a “ir) 


(n+2)(n + l)n (n—1) 
24 ‘ 


f I ' r (a)- 


Sumando miembro a miembro estas dos formulas, sacando 
factor comun a los terminos semejantes enambasi dividien- 
do por 2 todos sus terminos, obtendremos la siguiente for¬ 
mula: 













868 


MEMORIAS CIENTIFICAS I LITERARIAS 



(n + 2) (n + 1 ) n t n— 1 ) (n—2 ) /N (n + ‘2)(n + l) n> (n—1) (n—2) 

120 fv(aj + “720 
























Se puede notar en esta ultima formula que todas las diferencias pares tienen n elevado al cuadra 
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Pero para esplicar mejor el desarrollo deestas trasforma- 
ciones, daraos en seguida una eompleta para el coeficiente 
de la cuarta diferencia en la formula de Newton, la que pue- 
de servir de ejemplo para determinar de igual manera los 
otros valores. 

As i: 


n (n — 1) (n - 2) (n — 3) 0,12r — 0,88)( — 1,88) (— 2,88) 

24 — 24 


( —0,1056) ( —1,88) (— 2,88) fQ. 198528) (—2,88) 

24 ” 24 


— 0,57176064 
~24“ 


0,0238 


A la simple vista se puede notar que los coeficientes bino- 
mios de n son mucho mas pequeflos en la formula de Stir¬ 
ling que en las de Gauss.i Newton, lo que la hace mas apta 
que las anteriores, para los calculos de interpolacion, porque 
los errores introducidos seran tambien mas pequenos. 
Repitamos la formula de Gauss: 

f (a + nw) = f (a) + n f' | a +- 4- — ——— t” (a) 

+ ,4 a + - 1 -]+ 2 i tlv(a) 


liagamos 












DIFERENCIAS E INTERPOLAC10NES 


871 


La formula se trasformara en esta otra: 


a+ ■|') _f(a)+ ~2' f ( a+ ~2~) 


2 ' 2 
2 


f”(a) 


+ 


1 _ _ 1 _ 

2 ' 2 ' 2 


'• f ’” | a + "2”) + ’W f ‘ V(a) 


(A) 


La trasformacion del coeficiente binomio de la cuarta di- 
ferencia se hace de esta manera: 


1 L _ 1 _ 3 

<n+l)n(n—l)(n — 2) ~2 * ~2 ' ~h ' ~2 

24 ~ 8 



I haciendo la reduccion de los coeficientes binomios en 
los terminos de la formula (A) obtendremos la nueva for¬ 
mula: 


f a + 

\ 


w 


)-f(a) + -i r ( a +i)- f r 


(a) 


16 


' 2 


■]+ -128- f ' V(a) 


i haciendo 
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obtendremos: 


f ( a + "t) ~ f l a + w_ Tj“ f(a+1)_ T r ( a + t) 


-g- f” (a + 1) + 



a + 



3 

128 


f»v( a + l) 


Sumando miembro a miembro estas dos formulas i divi- 
diendo por 2 todos sus terminos obtendremos: 

,{ . w ) f(a) + ffa-r 1) 1 (f”(a) + f”(a-fl)) 

1 l a + ~r -2- 1 - 2 

3 f IV (a)+f IV (a+1) 

+ 128 2 


i sustituyendo respectivamente los factores: 


f(a) + f(a + l) . f” i 'a) + f”(a + l) _ fiv (a)+ fiv (a +1) 
2 ’ 2 y 2 


por sus correspondientes diferencias: 

resultara la^formula conocida con el nombre de «Formula de 
Stirling », cuya espresion jeneral es la que sigue: 
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Debemos hacer notar que en esta formula, la mas sencilla 
i comoda de todas para calcular i retener en la memoria, 
solo se emplean los promedios de las diferencias pares, que- 
dando por lo tanto todas las diferencias en una sola linea 
horizontal. 


Como complemento de estas formulas de interpolacion 
servira tambien el cuadro que sigue, el que nos dara una 
idea clara de la manera como se reparten los errores en las 
diferencias, cuando se ha cometido un error + x en la fun- 
cion orijinal f (a). 





Tabla de la reparticlon de I os errores 
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Se puede notar, ipso facto, que los errores se reparten, a 
partir de la fancion primitiva f (a), en forma de abanico, i 
que el mayor error se halla en la cuarta diferencia, en la 
misma linea de la funcion, lo que servira para indicar al 
momento' en cual de las funciones primitives se ha cometi- 
do una equivocation i proceder a correjirlo. 


EJERCICIOS FINALES 

P.— ENCONTRAR POR LAS FORMULAS DE INTERPOLACION DE 
NEWTON I DE GAUSS EL LOG. SEN DE 22° 24' 


Tabla de las diferencias 





Diferencias 



Argumento 

Funciones 






f 

r 


flY 

fv 



sen 20° 

9,53405 

+ 2028 





» 21° 

9,55433 

4 1925 

— 103 

■+8 



» 22° 

9,57358 

+ 1830 

— 95 

0 

0 

4 8 

» 23° 

9,59188 

+1743 

-87 

0 

—4 

+8 

» 24° 

9,60931 

— 79 

-4 



+1664 


+8 



» 25° 

9,62595 

+1589 

— 75 




» 26° 

9,64184 

1 






w=l° = 60' 
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log sen 22« 24’ «* 9,58101 
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log sen 22" 24’= 9,58101 
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Las tres formulas dan los mismos resultados. 
















